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Abstract 

We present a measurement of semi-inclusive spin asymmetries for positively and 
negatively charged hadrons from deep inelastic scattering of polarised muons on po- 
larised protons and deuterons in the range 0.003 < x < 0.7 and Q 2 >1 GeV 2 . Com- 
pared to our previous publication on this subject, with the new data the statistical 
errors have been reduced by nearly a factor of two. Prom these asymmetries and our 
inclusive spin asymmetries we determine the polarised quark distributions of valence 
quarks and non-strange sea quarks at Q 2 =W GeV 2 . The polarised u valence quark 
distribution, Au v (x), is positive and the polarisation increases with x. The polarised 
d valence quark distribution, Ad v (x), is negative and the non-strange sea distribu- 
tion, Aq(x), is consistent with zero over the measured range of x. We find for the 
first moments /„* Au v {x)dx = 0.77±0.10±0.08, /„* Ad v (x)dx = -0.52±0.14±0.09 
and Jq Aq(x)dx = 0.01 ± 0.04 ± 0.03, where we assumed Au(x) = Ad(x). We 
also determine for the first time the second moments of the valence distributions 
Jq 1 xAq v {x)dx. 
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Measurements of spin asymmetries in polarised deep inelastic scattering provide informa- 
tion about the spin structure of the nucleon. In particular detailed information can be 
obtained from polarised semi-inclusive deep inelastic scattering, jl+N —>■ fi + X + h, where 
in addition to the scattered lepton hadrons (h) are also detected. In this paper we present 
new results on semi-inclusive asymmetries from SMC data including those published in 
0. We analyse these asymmetries together with the inclusive asymmetries in the frame- 
work of the quark parton model (QPM) and determine the polarised quark distributions 
of the valence quarks and of the non-strange sea quarks. This is only possible due to the 
semi-inclusive data. 

Our experimental setup at the CERN muon beam consists of three major com- 
ponents: a polarised target, a magnetic spectrometer and a muon beam polarimeter. A 
detailed description of the experiment and the analysis of the inclusive data can be found 
in ref. 0, Positive muons of a nominal energy of 100 and 190 GeV were used. Most 
of the data were taken at 190 GeV. The muon beam polarisation Pb was determined 
from measured spin asymmetries in muon-electron scattering and for the 190 GeV data 
in addition from the energy spectrum of positrons from muon decays. It was found to be 
-0.795 ± 0.019 (-0.81 ± 0.03) for an average beam energy of 187.4 GeV (99.4 GeV). 
The target consisted of two cells filled with butanol, deuterated butanol or ammonia. 
The two cells were polarised in opposite directions by dynamic nuclear polarisation. The 
average polarisations Pt were approximately 0.90 for protons and 0.50 for deuterons with 
a relative error APt/Pt of 3%. The determination of semi-inclusive asymmetries requires 
the separation of hadrons from electrons which originate mainly from photon conversions. 
For this purpose we used a calorimeter [|J which consists of an electromagnetic and a 
hadronic part. The electromagnetic part amounts to 20 radiation lengths which is suf- 
ficient to contain electromagnetic showers. The total thickness of the calorimeter is 5.5 
nuclear interaction lengths. For each shower the ratio of the energy deposited in the elec- 
tromagnetic part to the total deposited energy is calculated. Electrons are eliminated by 
requiring this ratio to be smaller than 0.8. There was no hadron identification, only the 
charge of the hadron is known. 

The asymmetries of the spin-dependent virtual photon absorption cross sections for 
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production of positive (negative) hadrons are defined as 



(7 n -t- (7 TT 

where the indices || and ft refer to the relative orientation of the nucleus (proton or 
deuteron) and photon spins. Contributions from the asymmetry Ai are neglected. The 
cross sections cr + (~) refer to the number of produced hadrons, i.e. all hadrons detected in 
one event are counted. The extraction of the semi-inclusive asymmetries from the counting 
rates is described in |L|. 

The counting rate asymmetry A exp is related to the virtual photon asymmetry A\ 

by: 

The depolarisation factor D depends on the event kinematics and on the ratio R of 
longitudinal to transverse virtual photon cross sections 0] . The dilution factor / accounts 
for the presence of unpolarisable nuclei in addition to the protons or deuterons in the 
target. It can be expressed in terms of the numbers ua of nuclei with mass number A and 
the corresponding total spin-independent cross sections a 1 ^ ■ Taking into account radiative 
effects on the nucleon [pi the effective dilution factor is defined as 



I7 tot 

r _ a p,d n p,d a p,d /q\ 
J ~ n tot y- tot ■ \ 6 ) 



Here p, d stands for proton or deuteron. 

The evaluation of the dilution factor for inclusive events, /*, is described in ref. J2]. 
The main contributions to the total inclusive cross section are the one photon exchange 
with vertex and vacuum polarisation corrections, and the inelastic, elastic and quasielas- 
tic bremsstrahlung processes. In the elastic and quasielastic bremsstrahlung processes no 
hadrons are produced and therefore these processes do not contribute to the total semi- 
inclusive cross section. In the analysis hadrons are selected if their energy is above a 
certain fraction, z, of energy transfer v. This reduces the energy available for a possible 
bremsstrahlung photon accompanying an inelastic event, relative to the inclusive case. 
The contribution from inelastic bremsstrahlung is then reduced. This reduction was es- 
timated to be small and is included only in the systematic error of the dilution factor. 
We compute the dilution factor for the semi-inclusive events, f sl , without the elastic and 
quasielastic contributions in the cross section. In the inclusive case, at low x there is a 
large contribution of elastic bremsstrahlung from high Z nuclei in the denominator of 
eq. ([I]). It is not present in the semi-inclusive case and therefore f st is 30% larger than f l 
in this region, whereas the difference is very small for x > 0.1. Using the QPM we com- 
pute from f st the dilution factors for positive, / + , and negative, /~, hadrons || 0. For 
the polarised deuterated butanol target / + and /~ are equal to f st . For the butanol and 
ammonia targets, where protons are polarised, the ratio f + / f st ~ 1.07 and f~ / f s% ~ 0.88 
at x = 0.5 and both ratios are close to 1 at low x. The ratios differ from unity for these 
targets because more positive hadrons are produced on protons at large x as compared 
to the isoscalar unpolarised target nuclei. The correctness of the above procedure for 
the evaluation of the dilution factor was verified using a Monte Carlo in which radiative 
processes and hadron production were simulated. 
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Polarised radiative corrections are applied to the asymmetries as described in ref. |2], 
[8| . In this procedure, as for the evaluation of the dilution factor, contributions of processes 
where no hadron is produced are omitted. 

In order to interpret the results in the QPM a cut of Q 2 > 1 GeV 2 is applied. The 
current fragmentation region is selected by a cut oia z = Ehad/v > 0.2. After cuts on u, y, 
the energy and the angle of the scattered muon as in ref. 0, we obtain 32 x 10 6 events. 
After applying the z cut we obtain 5 x 10 6 positive and 4 x 10 6 negative hadrons. All of 
the data are in the deep inelastic region, W > 3 GeV, and cover the range 0.003 < x < 0.7 
with an average Q 2 ~ 10 GeV 2 . 

The asymmetries measured in different periods of data taking are compatible with 
each other, and thus all SMC data are combined. The semi-inclusive asymmetries ~ d for 
positive and negative hadrons from the deuteron and the proton are presented in Fig. |] 
and Tab. [I]. The correlations between the asymmetries are listed in Tab. The main 
contributions to systematic errors are due to the uncertainties of the target and beam 
polarisations, to the variation in time of the spectrometer acceptance and to the uncer- 
tainty of the dilution factor. Negligible contributions arise from secondary interactions of 
hadrons in the target, radiative corrections and differences in acceptance for pions, kaons 
and protons due to different absorptions in the target and different angular distributions. 

The inclusive asymmetry A\ v can be expressed in terms of polarised and unpolarised 
structure functions: 

A ^-w^rwB) 2x ^ +R{x ' Q ^- (4) 

In the QPM the structure functions can be written in terms of polarised quark distri- 
butions Aq = — q* and unpolarised quark distributions q = q^ + q^, where q*(q^) 
is the distribution of quarks of flavour q and spin parallel (antiparallel) to the nucleon 
spin. The published leading order parametrisations of the unpolarised quark distribu- 
tions were obtained from experimental values of F$ using the relation F%{x, Q 2 ) = 
x J2 q e l [q( x , Q 2 ) + Q~( x , Q 2 )] ■ The values of Ff were extracted from cross sections assuming 
non-zero values of R. Therefore the expression for Ai p becomes: 

lp( *' g) " E q e 2 [q(x,Q 2 ) + q(x,Q 2 )] L 1 ^ 

where e q is the fractional charge of the quark q with q = u,d,s. For x < 0.12 we use 
a parametrisation of R measured by the NMC [[]] and for higher x we use the SLAC 
parametrisation At our average Q 2 of 10 GeV 2 the ratio R varies from 0.15 at 

x = 0.005 to 0.07 at x = 0.5. 

For semi-inclusive processes the asymmetries depend in addition on the fragmenta- 
tion functions D^(z, Q 2 ), which represent the probability that a quark q fragments into a 
hadron h. Using the quantity D^(Q 2 ) defined as D^(Q 2 ) = J^ 2 dz D 1 ^ (z, Q 2 ) we have: 

A H-)r n ^ [Ag(x,Q 2 )DZ(Q 2 ) + ^,Q 2 )Dl(Q 2 )} r } 

A^ '(x,Q) = L — r l + R{x,Q*) , (6) 

E q , h e 2 q [q(x,Q 2 )D^) + q(x,Q 2 )Dl(Q^} [ ^ 

where the summation is over 7r + , K + and p for positive hadrons and over 7r~, K~ and p 
for negative hadrons. We assumed the same R for positive and negative hadrons, and used 
the same R as for the inclusive case. It was checked that our acceptance depends only 
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weakly on z in the region covered by the integral of the fragmentation functions. Since the 
fragmentation process is parity conserving and the polarisation of hadrons is not observed, 
the fragmentation functions do not depend on the quark helicity (D q i = D q i). In general 

D^, so the measurement of semi-inclusive asymmetries allows a separation of Aq and 
Aq, whereas only the sum, (Aq + Aq), can be determined from the inclusive asymmetries. 

Using isospin symmetry, similar expressions are obtained for the deuteron asymme- 
tries. From the data on proton and deuteron we can separate spin contributions from u 
and d quarks. The spin dependent deuteron cross section is considered to be the sum of 
the proton and the neutron cross sections with a correction to account for the .D-state 
probability ojd = 0.05 ± 0.01 of the deuteron, as in ref. ||. We use the same R for the 
proton and the deuteron. 

With eqs. and (^) and the corresponding relations for the deuteron our measured 
inclusive and semi-inclusive asymmetries can be used to evaluate the polarised quark 
distributions. The polarised valence quark distributions are defined as Au v (x) = Au(x) — 
Au(x) and Ad v (x) = Ad(x) — Ad(x). To reduce the number of unknowns we assume a 
SU(3)f symmetric sea which we denote by Aq(x): 

Aq(x) = Au(x) = Ad(x) = As(x) = As(x) . (7) 

The sensitivity of our results to this assumption will be discussed below. In particular the 
assumption involving As(x) and As(x) has a negligible influence on our results, because 
our hadron sample contains mainly pions, thus Aq(x) provides information mainly about 
non strange sea quarks. 

The six measured asymmetries are linear combinations of three unknowns: Au v , 
Ad v and Aq. Six equations for the asymmetries can thus be written in matrix form: 

A = BAq (8) 

where A = (A lp , Af p , A^ p , Ai d , A^ d , A~[^j and Aq = {Au v , Ad v , Acfj. The elements of the 
matrix B are determined from eqs. (|5D and (H). They depend on the unpolarised quark 
distributions q, the fragmentation functions D^, the ratio R, and up. 

The asymmetries were measured at average Q 2 values in each x bin varying from 
1.3 GeV 2 at our lowest x to 60 GeV 2 at x = 0.5. No significant Q 2 dependence is observed 
in the inclusive data 0, [11], therefore we assume the asymmetries at measured Q 2 to be 
equal to those at Q 2 = 10 GeV 2 . In eq. (H) we use parametrisations of the unpolarised 
quark distributions |I2 |, of the ratio R, and of the fragmentation functions all at 



Q 2 = 10 GeV 2 . For the unmeasured fragmentation functions additional assumptions are 
made, as discussed below. 

In order to determine the leading order polarised quark distributions Au v (x), Ad v (x) 
and Aq(x) we solve eq. (|]) in every x-bin independently by minimising the quantity 

X 2 = (A — BAq) T (C0V4)- 1 (A- BAq) , (9) 

where Cov^ is the covariance matrix of the asymmetries (cf. Tab. [|). The results are shown 
in Fig. ^| and Tab. We see in the figure that the statistical error on Aq(x) in the region 
x > 0.3 is larger than q(x) which is an upper limit for |Ag(x)|. In this region, in order to 
reduce the statistical error on Au v (x) and Ad v (x), we set Aq(x) = before solving the 
system of equations. The results are shown as open circles in Fig. 2. Uncertainties due to 
this assumption were included in the systematic error. 
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We observe that Au v (x) is positive while Ad v (x) is negative. The polarised sea 
quark distribution Aq(x) is compatible with zero over the full range of x. Fig. |3] shows the 
polarisation Aq v (x) / q v (x) of the valence quarks. The average polarisation is approximately 
50% for the u v quarks and —50% for the d v quarks. For the u v quarks we observe a lower 
average polarisation of 0.18 ± 0.10 in the interval 0.003 < x < 0.06 (first six bins) and a 
higher average polarisation of 0.57 ± 0.05 in the interval 0.06 < x < 0.7 (last six bins). 

The uncertainties of these results are dominated by the statistical errors. Contri- 
butions from different sources to the systematic errors on the integrals / ° 7 03 Aq(x)dx are 
listed in Tab. f| and the values of the integrals in Tab. |]. The largest contribution is re- 
lated to the assumptions made for the unmeasured fragmentation functions of s quarks. 
Strange quark fragmentation functions to pions were assumed to be equal to the measured 
unfavoured fragmentation function: Dg^ ,7T = . For the favoured fragmentation func- 
tions of s quarks to kaons we assumed Df = = . This is motivated by the fact 
that in both fragmentation processes, s —>■ K~ and d —>■ n~ , a im-pair has to be produced 
to create the hadron. Values of the unmeasured fragmentation functions obtained with 
the above assumptions were varied by a factor 0.5 to 1.5 and the difference in the resulting 
polarised quark distributions was included in the systematic error. The uncertainty due 
to the unpolarised quark distributions was taken to be the difference resulting from two 
different parametrisations, GRV-94 (LO) @ and CTEQ 3L (LO) 0. The systematic 
error due to the assumption Aq(x) = in the last two x bins was calculated by setting 
Aq(x) = ±q(x) in these bins. 

In the QPM, with the assumption Au(x) = Ad(x), the quantities 6[g{(x) — <7i(x)] 
and Au v (x) — Ad v (x) are equal. The former is obtained from inclusive data only while the 
latter can be extracted using only the four semi- inclusive asymmetries in eq. (H). The two 
quantities are generally in good agreement, as seen in Fig.[|, showing the consistency of 
the present analysis. The discrepancy in the last bin is mainly due to the low value of the 
asymmetry Af in this bin. This is also reflected in Tab. |3] by the low x 2 probability for 
this bin. It has been checked that this low value could not be caused by a false asymmetry 
due to a change in the spectrometer acceptance. With the constraint Aq(x) = only two 
unknowns, Au v (x) and Ad v (x), have to be determined. This can be done using only the 
inclusive asymmetries A\ v and Aid which have a better statistical accuracy. For this reason 
adding the semi-inclusive asymmetries has very little impact in the last two bins, i.e. the 
results for Au v (x) and Ad v (x) from A\ p and Aid only are almost identical with the values 
shown in Tab. [3] where all six asymmetries are used. 

In order to calculate first moments, Aq = f Aq(x)dx, of the polarised quark distri- 
butions, we need extrapolations of Aq(x) to the unmeasured regions of x. Contributions 
from the unmeasured large x region (0.7 < x < 1) are negligible because their upper 
limits given by the unpolarised distributions are small there. The low x (0 < x < 0.003) 
extrapolation was done in several ways. Two functional forms of the polarised quark dis- 
tributions were fitted to our semi-inclusive and inclusive asymmetries J7j. We also used 
parametrisations of the polarised quark distributions obtained from two different QCD 
analyses |15], [l6j of inclusive data. For every quark distribution we obtained four values 



for the integral / 003 Aq(x) and estimated the mean value and the error. For the valence 
quarks these values (see Tab. |5|) are much smaller than the upper limits given by the un- 
polarised distributions | J a003 Aq v (x)dx\ < J 003 q v (x)dx. This upper limit is 0.22(0.17) 
for GRV-94, LO parametrisation of the u v (d v ) quarks. There is no such constraint for the 
sea quarks. The first moments of the polarised quark distributions are 
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Au v = 0.77 ±0.10 ±0.08, 
Ad v = -0.52 ±0.14 ±0.09, 
Aq = 0.01 ±0.04 ±0.03. 

The statistical errors are reduced by almost a factor of two compared to our previous 
publication [lj . The contributions to the moments from measured and unmeasured regions 
are detailed in Tab. ||. 

In order to study the effect of the assumption of a SU (3) f symmetric sea in eq. (|7|) 
we replace it by only an isospin symmetric sea, Au(x) = Ad(x) = Aq(x) and allow 
As(x) = As(x) = i]Aq(x), varying rj in the range 0.25 < r\ < 1.5. The corresponding 
variations of the first moments are of the order of 0.01. Similar variations are observed 
if we use As(x) as given in ref. jlS] , where the integral is Jq 1 As(x)dx = —0.06, a value 
consistent with the result Jq 1 (As(x) + As(x))dx = — 0.12±0.04 obtained from the inclusive 
analysis [JO ] . This shows that the measurement is only sensitive to the polarisation of the 
non-strange sea quarks. 

Releasing the constraint Au(x) = Ad(x) and setting As(x) = As(x) = (Au(x) + 
Ad(x))/2 leads to the results given in the lower part of Tab. |. The results are consis- 
tent with those obtained with the constraint Am = Ad. The statistical error on Au v (x) 
practically does not change, whereas the error on Ad v (x) increases by a factor of two. 

Polarised semi-inclusive deep inelastic scattering is presently a unique tool to mea- 
sure the polarised valence and sea quark distributions. However, the corresponding first 
moments can be determined by other methods. With the assumption of a SU (3) / symmet- 
ric sea (eq. (|7|)) and SU(3)f symmetry for the weak decays in the baryon octet, the first 
moments of the valence distributions can be obtained from the axial matrix elements, F 
and D, of baryon decays. It should be noted that the results depend strongly on these two 
assumptions, whereas the results from our semi-inclusive analysis are insensitive to any 
SU(3)f assumptions, as discussed above. With the values quoted in Ref. |L7|, |18| one finds 



Au v = 2F = 0.92 ± 0.02 and Ad v = F - D = -0.34 ± 0.02, which is in good agreement 
with our results of Tab. |5| obtained from deep inelastic scattering only. Our result for Aq 
is consistent with the first moment of the polarised strange quark distribution obtained 
from the inclusive analysis JT1J quoted above. 



From our polarised quark distributions we find for the second moments: 

f 1 xAu v (x)dx = 0.155 ± 0.017 ± 0.010 and C xAd v (x)dx = -0.056 ± 0.026 ± 0.011 , 
Jo Jo 

where the first error is statistical and the second is systematic. Due to the additional factor 
x the low x extrapolations to the unmeasured region give negligible contributions while 
the large x contributions remain small. A calculation of the second moments in lattice 
QCD @ predicts /* xAu v (x)dx = 0.189 ± 0.008 and J* xAd v (x)dx = -0.0455 ± 0.0032, 
which is in good agreement with our measured values. 

To summarise, we have measured semi-inclusive spin asymmetries for positively and 
negatively charged hadrons from polarised protons and deuterons and we have determined 
the polarised quark distributions Au v (x), Ad v (x) and Aq(x) of the valence and the non- 
strange sea quarks. The polarisation of the u v quarks was found to be positive and to 
increase with x. The polarisation of the d v quarks was found to be negative while the 
polarisation of the non-strange sea quarks is consistent with zero. Our results for the first 
moments of the polarised valence quark distributions are in good agreement with results 
obtained from the axial matrix elements F and D. The second moments are consistent 
with lattice QCD predictions. 
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Table 1: Values of semi-inclusive spin asymmetries Af p , A lp , A\ d and A ld . The first error 
is statistical, the second one systematic. 
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Table 2: The correlation coefficients between the asymmetries, tn the correla- 
tions only the dominating statistical uncertainties are taken into account. There- 
fore the proton and deuteron asymmetries are uncorrelated. The correlations are 
given by cor(A+, A~) =< n + n~ > /V< n+ 2 >< rr' 1 > and cor(A, A +( -A) =< > / 

\J < n+( - ) 2 >, where n +( -^ is the number of positive (negative) hadrons per scattered 
muon. 
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Table 3: Values of the polarised quark distributions xAu v (x), xAd v (x) and xAq(x). The 
distributions were obtained with the assumption Au(x) = Ad(x). The first errors are 
statistical and the second ones are systematic. The values in the last two bins of x cor- 
respond to the open circles in Fig. [J. In the last two columns the value of x 2 / n an d the 
corresponding probability p are given. 
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Table 4: Contributions to the error on the integral J 003 Aq(x)dx. The abbreviation f.f. 
stands for fragmentation functions. The values of the integrals are given in Tab.[|. 
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Contributions to the first moments 




of the polarised quark distributions 


Au(x) = Ad(x) 


X 


0-0.003 


0.003 - 0.7 


0-1 


Au v 


0.04 ±0.04 


0.73 ±0.10 ±0.07 


0.77 ±0.10 ±0.08 


Ad v 


-0.05 ±0.05 


-0.47 ±0.14 ±0.08 


-0.52 ±0.14 ±0.09 


X 


0-0.003 


0.003 - 0.3 


0-1 


Aq 


0.00 ±0.02 


0.01 ±0.04 ±0.03 


0.01 ±0.04 ±0.03 


Au(x) ^ Ad(x) 


X 


0-0.003 


0.003 - 0.7 


0-1 


Au v 


0.04 ±0.04 


0.72 ±0.11 ±0.06 


0.76 ±0.11 ±0.07 


Ad v 


-0.05 ±0.05 


-0.45 ±0.30 ±0.25 


-0.50 ±0.30 ±0.25 


X 


0-0.003 


0.003 - 0.3 


0-1 


Au 


0.00 ±0.02 


0.01 ±0.05 ±0.01 


0.01 ±0.05 ±0.02 


Ad 


0.00 ±0.02 


0.01 ±0.14 ±0.12 


0.01 ±0.14 ±0.12 



Table 5: The values for the extrapolations to low x (first column), the integrals over the 
measured range (second column) and the first moments (third column). In the last two 
columns the first error is statistical, the second systematic. The upper part of the table 
shows the results obtained assuming Au(x) = Ad(x), the lower part shows the results 
without this assumption. 
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Figure 1: Semi-inclusive spin asymmetries for the proton and the deuteron as a function 
of x for positively and negatively charged hadrons. The error bars are statistical and the 
shaded areas represent the systematic uncertainties. 
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Figure 2: The polarised quark distributions xAu v (x), xAd v (x) and xAq(x) obtained with 
the assumption Au(x) = Ad(x). The open circles are obtained when the sea polarisation 
is set to zero while the closed circles are obtained without this assumption. The error 
bars are statistical and the shaded areas represent the systematic uncertainty. The curves 
correspond to the upper and the lower limits ±xq(x) from the unpolarised quark distribu- 
tions evaluated at Q 2 = 10 GeV 2 . In the bottom plot the curves are ±x(u(x)+d(x))/2. 
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Figure 3: The polarisation Aq v (x) / q v (x) of the valence quarks. Only the statistical errors 
are given. 
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Figure 4: Tie difference between the spin-dependent structure functions of the proton 
and the neutron, 6x[gl(x) — g™{x)\, (open circles), and the difference between polarised 
up and down valence quark distributions x[Au v (x) — Ad v (x)], as determined with the 
assumption Au(x) = Ad(x) from SMC semi-inclusive asymmetries only (closed circles). 
Errors are statistical. 
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